Strength Distribution in Gradient Networks 
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This article describes a gradient complex network model whose weights are proportional to the 
diflFerence between uniformly distributed "fitness" values assigned to the nodes. It is shown analyti- 
cally and experimentally that the strength (i.e. the weighted node degree) density of such a network 
model can be well approximated by a power law with 7 « 0.35. Possible implications for neuronal 
networks topology and dynamics are also discussed. 
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I. INTRODUCTION 

Great part of the interest focused on complex networks 
Q 13 stems from scale free or power law distributions 
of respective topological measurements, such as the node 
degree. At the same time, weighted complex networks 
have attracted growing interest because of their relevance 
as models of several natural phenomena, with special at- 
tention to systems used for distribution/collection of ma- 
terials or information. 

There are two main ways to approach the degrees of 
weighted networks: (i) by thresholding the weights and 
using the traditional node degree and (ii) by adding 
the weights of the edges attached to each node, yielding 
the respective node strength Related recent works 
include the identification of strength power law in word 
association networks |^ 1^ , studies of scientific collabora- 
tions and air-transportation networks (3] , the analysis of 
amino acid sequences in terms of weighted networks 6] , 
the investigation of weighted networks defined by dynam- 
ical coupling between topology and weights |4|L analyti- 
cal characterization of thresholded networks as well 
as the characterization of motifs and shortest paths 
in weighted networks 10]. 

The current article describes a gradient oriented net- 
work whose nodes have a respective "fitness" value vary- 
ing uniformly between and 1 with resolution A^, and 
every node is connected to all other nodes with higher 
fitness through an edge with weight proportional to the 
difference between the nodes fitness. It is shown both an- 
alytically and experimentally that the strength densities 
obtained for gradient networks can be well approximated 
by a power law with parameter 7 « 0.35. 

We start by describing the gradient network model and 
follows by calculating its respective strength distribution, 
illustrating experimentally obtained dilog curves, and 
discussing possible implications for neuronal networks. 



II. GRADIENT NETWORKS 

Let the initial state of the directed network F contain 
N nodes such that each node i is potentially connected 



to every other node j, therefore implying a maximum 
of N"^ - N edges. A "fitness" value v{i) = [aNr\/Nr — 
where a is uniformly distributed in the interval [0, 1) and 
Nr is a positive integer parameter expressing the total of 
possible different values v{i) — is randomly assigned to 
each node i of the network. Consequently, the values 
v{i) are equally spaced, i.e. v{i) — qAr for some q S 
{0,l,...,Nr- 1}, Ar = 1/Nr and < v{i) < 1. For 
N > Nr, the expected number of nodes with each specific 
value v{i) can be estimated as n = NAr. The weight of 
the edge connecting node i to node j, namely w{i,j), with 
i 7^ j, is defined as in Equation It follows that < 
w{i,j) < 1. Each weight u;(i, j) is henceforth represented 
as the entry W{j, i) of the corresponding weight matrix 
W. 



w{i,j) = 



{v{i) - v{j))Ar ifv{i) > v{j) 







otherwise 



(1) 



Figure^illustrates a simple gradient network obtained 
by the above described procedure. Observe that the 
weights can be understood as being proportional to the 
gradient of v{i) while moving from any of the neighbors 
j of i, such that v{j) < v{i), to node i, hence the name 
gradient networks. 



III. ANALYTICAL DERIVATION OF THE 
INSTRENGTH DENSITIES 



The node instrength of a node i belonging to a directed 
weighted network with weight matrix W is defined as 
corresponding to the sum of the weights of the incoming 
edges, as expressed in Equation [3 Considering a node i 
with specific "fitness" v{i) = qAr, we expect to find, in 
the average, nq network nodes with fitness smaller than 
v{i), namely the nodes having "fitness" v{i) — jAr with 
j = 0, 1, . . . , (7 — 1 and respective weights {q — j)A'^ — 
pA^, where p ~ q — j — 1,2, ... ,q. The instrength of 
node i such that v{i) — qAr can thus be estimated as in 
Equation 121 
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FIG. 1: A simple network obtained by applying the proce- 
dure described in the text. Observe that each node receives 
connections from nodes with smaller fitness values, while the 
respective weights are proportional to the fitness diflterences. 



By series expansion around s = it follows that 



log{h{s))Ki a--fs + (3s^ + 0{s^) (8) 

where 7 = -0.25(V^ - 2)/(x/2 - 1) w 0.354, (3 = 
1/32(4- 3%/2)/(%/2- 1)2 -0.044 and a is a function 
of N, A and A^. Since the coefficient of the term in 
is about 10 times smaller than the terms in s, the dilog 
curve of the instrength distribution estimated from the 
respective histogram for A 3> A^ can be well approxi- 
mated by a straight line with slope —7. 



N 

%)=^W^(*,j) (2) 
i=i 

k{q) w n 2^ pA^ = (3) 

p=l 

By making c = 0.5nA^ — 0.5A^A^, we have 

fc(g) « ciq^ + q) (4) 

The maximum degree is therefore obtained for v{i) = 
{Nr-l)Ar, yield ing k^nax ~ 0.5y2. The fact that we know 
from the network construction that a different instrength 
value is obtained for each integer instance q, allied to 
the fact that in the average a total of n nodes will have 
v(i) = qAr, imply that the continuous density (over the 
free variable x) of instrengths is given as in Equation El 

p{x) = 6{x~ k{q)) (5) 

Let us now consider the estimation of the indstrength 
distribution by using the traditional binning scheme 
given by Equation which maps each indegree k into 
the respective bin of length A indexed by K. Observe 
that symmetric rounding can also be considered by ex- 
changing the floor function by the round function. 

K = [fc/AJ (6) 

Such a binning procedure can be understood as parti- 
tioning the X-axis with interval size A. By approximating 
the values of k given in Equation 01 as fc w cq^ , we have 
that q = \fkjc and the number of cases expected in each 
bin — yielding the histogram h{K) — can be estimated 
as 

h{K) « n (V(i^ + l)A/c - y/KKf^ (7) 



IV. DISCUSSION 

Figure [3 shows the dilog plot of instrength density ob- 
tained after 1000 simulations of a gradient network as- 
suming N = 2000, Dr = 0.001 and A = 0.01. As ex- 
pected, the resulting curve closely resembles a straight 
line with estimated (linear regression) absolute slope 
value equal to 0.41 « 7. 

We have described a simple gradient network model 
and shown that its node instrength distribution is close 
to a power law with parameter 7 « 0.35. This result im- 
plies that although several nodes of a gradient network 
have low strength values, there are hubs characterized by 
varying large instrength values. Because the instrength of 
a node can be understood as the total weighted influence 
it receives from the adjacent nodes, it follows that in case 
the node state is directly proportional to its instrength, 
the dynamics of such networks will be characterized by 
a near power law with similar parameters as those of the 
instrength density. Such an interpretation is particularly 
interesting from the perspective of integrating neuronal 
network and complex network research (e.g. |lll |. where 
each neuron is represented by a node and the synaptic 
connections by directed edges whose weights reflect the 
respective synaptic strengths. The "fltness" values in 
such a case could be related to gradient of neurotrophic 
growth factors or depolarization bias facilitating the ac- 
tion potential [I3, • 

It would be interesting to investigate the strength dis- 
tribution for "fltness" statistical models other than that 
presently considered uniform density. 
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FIG. 2: Dilog diagram (average ± standard deviation of instrength density (non-cumulative) estimated from the respective 
histogram by using N = 2000, Ar = 0.001 and A = 0.01. An almost straight line segment curve is obtained with slope —0.41. 
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